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ABSTRACT 

Using a generally accepted model we present a comprehensive analysis (within the page limitation) of an Euler- 
Bemoulli beam with PZT sensor-actuator and pure rate feedback. The emphasis is on the root locus — the 
dependence of the attainable damping on the feedback gain. There is a critical value of the gain beyond which the 
damping decreases to zero. We construct the time-domain response using semigroup theory, and show that the 
eigenfunctions form a Riesz basis, leading to a “modal” expansion. 


1. INTRODUCTION 

In this paper we present a comprehensive analysis of an Euler-Bemoulli beam with PZT sensor-actuator along its 
entire length. The sensor output is a charge in a condenser and the actuator input is the current, a differentiator 
circuit being then an essential component, yielding “rate feedback.” We use a generally accepted model. Tzou 
et a/. 4 present purely computational results and seem to be unaware of a purely theoretical analysis given earlier by 
Chen et al A The most important design parameter is the control gain and the damping attainable — we construct a 
full root-locus analysis (omitting details to keep within the page limitation). We also unearth a curious phenomenon 
— the existence of a deadbeat mode (real eigenvalue) not noticed hitherto. We show that the eigenvalues are the 
roots of an entire function of order one-half, proving in particular the existence of a countably infinite number of 
eigenvalues. We also show that the eigenfunctions form a Riesz basis. We also construct the Green s function for 
the nonhomogeneous eigenvalue problem. As in Chen et al. we use the theory of semigroups of operators to obtain 
the time-domain solution. Our proof of the exponential stability is different from that in Chen et a/., as is our 
choice of the function space. We note that a similar analysis for a Timoshenko model (a “smart string” ) is given 
in Balakrishnan, 6 where there is a critical value of the gain at which there are no eigenvalues and the semigroup is 
actually nilpotent (“disappearing” solution). 


2. MAIN RESULTS 


The Euler-Bemoulli model formulates as 

cf""(t,s) + mf\t l 3) = 0, 0< s< L, 0< t 'I 

/(t, 0) = 0 = /'(t, 0); f"'(t,L) = 0 > (1) 

cf"(t,L)+af'{t l L) = 0 J 

where f(t,s) is the displacement and the superdots indicate derivative with respect to t and the primes indicate 
derivative with respect to s. It is convenient to set 


For a precise formulation of the time-domain response we need to specify first the choice of function spaces. We 
pick L 2 (0, L] for /(£, •). Let A a denote the operator defined by 

A 0 f = of"" 


where 


V(A 0 ) = [/ | € L 2 [0, L); /( 0) = 0 = /'(0) = f"'(L)} . 
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Let 


n = l 2 [o,l\ x e 1 . 

Define the operator A with domain and range in H by: 


Ax = 


A 0 f 

cf"(L) 


with domain 

V(A) = 

It is convenient to adopt the notation 


/ € V{A 0 ) and b= f'{L) 
At,x = cf"(L), x € V(A). 


Then for x in T>(A): 


[Ax,x] = f cf""(s)f( 

Jo 

= c[ L \f”( 

Jo 


s) ds + cf"(L) f'(L) 


s)\ 2 ds. 

It is readily seen that A has dense domain and is self-adjoint and nonnegative definite, and has compact resolvent. 
Also zero is not an eigenvalue. Let s/~A denote the positive square root. On the product space 


V (VI) x L 2 [0,L| 


introduce the “energy” inner product 


\Y,Z] e = [v/Iyi, v/Izi] + m (j/2, 22] 


Y = 


V1 , = 

*1 

2/2 

*2 


P(7I) 

For 1/1 in V(A ), we see that 


f"€L 2 [0,L} and b=f'(L), /(0) = /'(0) = 0 


[Y,Y} b = [A 2 / 1 , 2 / 1 ] + rn [ 2 / 2 . 2 / 2 ] ~ “energy” (potential + kinetic). 

We denote the product space under this inner product by He and note that it is a Hilbert space. Let A denote 
the operator defined by: 



/2 




/l() 

AY = 

-cK’(L) 

a 

Y = 

X 

h 

= 

/ m 


- Aoh 
m 



M-) 


and 

V{A) = 

Thus defined we can verify that 


/i(-) 

f[(L) 


€V(A) 

V(A) = T>(A') 


M-) 

fl'(L) 


v(Va). 


and that A is dissipative: 


h[{A + A*)Y, T] = Re[Ay,y]£: = ^||A,,x|| 2 = ^ c 2 |/i'(L)| 2 . 
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It is readily verified that A has a compact resolvent and that A generates a C 0 contraction semigroup. With these 
definitions, the system (1) goes over into the abstract formulation. 

V'(t) = •AV'(t). 

This choice of the function space is technically different from that in Chen et al . 5 

Eigenvalues and eigenfunctions of A . 

Our primary interest is in the modal decomposition — the eigenvalues of A and the corresponding eigenfunctions. 

Or, equivalently, in the resolvent of A. Let ft(A, .4) denote the resolvent of A. Let 

Tl(\,A)Y = Z 


where 


Since Z € T>{A), we can write 


and 


Y = 


*1 

b 

h 2 


Z 


/i(') 

f'i(L) 

/ 2 () 


(A I -A)Z = Y 


yields 


Hence 


A/ i - h = h i 

Apfi 
m 

cfW _ 


A/2 + 


- 


A/|(L) + 


= b. 


a 


av/im + /™ 

Au/i(L) + cf'i(L) 

MO) = o 


i/ 2 (/i 2 (s) + AM*)). 

ab 

/i(0) = 


0 < s < L \ 


( 3 ) 


Eigenvalues 

First we consider the eigenvalue problem, setting 

h\ — 0 = h,2\ 6=0. 


Let 




7 4 = -A V 


where 


A = |A|e ,e . 


Then the solution satisfying the conditions at zero yields: 


/,(«) = a(Cosh 7 s - Cos 7 s) + b(Sinh 7 s - Sin 7 s), 0 < s < L. 


The constants a and b are then determined by the conditions at L: 

a(Aa 7 (Sinh 7 L + Sin 7 L) + c 7 2 (Cosh 7 L + Cos 7 L)) + b(Aa 7 (Cosh 7 L - Cos 7 L) + c 7 2 (Sinh 7 L + Sm 7 L)) 

a 7 3 (Sinh 7 L — Sin 7 L) + b 7 3 (Cosh 7 L + Cos 7 L) = 0. 


0 
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Let 


H(\) = 


Aa7(Sinh7 L 4 - Sin7L) 

4 c7 2 (Cosh 7L 4 Cos 7L ) 


Aa7(Cosh 7 L — Cos 7L) 

4 C7 2 (Sinh 7L 4 Sin 7L) 


and 

Then 


7 3 (Sinh7L — Sin7L) 7 3 (Cosh7L 4 C0S7L) 


D( A) = Det H (A). 


( 4 ) 


D( A) = (7 4 ) [(Cosh 7L 4 Cos 7L)(Aa(Sinh 7L 4 Sin 7L) 4 c7(Cosh 7L 4 Cos 7L)) 

- (Sinh7L - Sin7L)(Aa(Cosh7L — C0S7L) 4 C7(Sinh 7L 4 Sin 7^))] 

= 27 4 [07(1 4 Cosh7L C0S7L) 4 Aa(Sinh 7L Cos 7L 4 Cosh7L Sin7L)] . ( 5 ) 

We note that zero is not an eigenvalue. The eigenvalues {A*} are thus determined by the nonzero roots of 
07(1 4 Cosh7L C0S7L) 4 Aa (Sinh7L Cos 7L 4 Cosh7L Sin7L) = 0 


Or, using 


we have 


A = 


-iT 


(1 + Cosh 7L Cos 7L) - — (Sinh 7L Cos 7L + Cosh7L Sin 7L) = 0 . 

vc 


( 6 ) 


Theorem 2.1^ 

A has exactly one real-valued eigenvalue. 


Proof 

Setting L ~ 1 , and using a to denote £ , and expressing the trigonometric products in (6) as sums, we have 

/ = 1 4 Cosh 7 Cos 7 — iQ7(Sinh7 Cos 7 4 Cosh 7 Sin 7) 

/ = 1 4 ^(Cos7(l 4 2) 4 Cos7(l - i)) 

- iQ7^[Sinh7(l + i) + Sin7(l + i) + Sinh7(l-») + Sin 7(1 — i)] . 

Hence making the 1:1 transformation 

7 = x(i — 1) 

we obtain 

/(^) = p(x) = 14 ^(Cos 2 x 4 Cosh 2 x) — ax(Sin 2 x 4 Sinh 2 x) ( 7 ) 

yielding an equivalent expression for determining the eigenvalues. Note that ^(*) is real-valued for real values of x. 
Further 

g(0) = 2 

while, as x — ► oo, (x real), we note that 

g(x) — ► oo. 

Hence there is a positive real root. Denote it xo- Then 

A = -7 2 i = —x 2 (i — 1 ) 2 i = — 2x 2 . 


Hence 

A 0 = - 2 xq 


1“Due to J. Lin; private communication. 
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is an eigenvalue. We note that x\ is the ony 
A , we must have, denoting it by A x , 


real-valued root of g (- ). Indeed, if there is a real-valued eigenvalue of 
A i = — 2x i 


and x 0 must be a root of p(-). Hence 


£Q = x l- 


Or, Ao is the only real-valued eigenvalue of A. 

We note that the corresponding eigenfunction is given by 

4, (a) = (Cosh 7 0 - Cos 7o ) (Cosh 7o s - Cos 7o s) - (Sinh 7o - Sin 7o )(Sinh 7o s - Sin 7o a) 


where 

7o = x o{i — 1)» A = — 2x 0 . 


Theorem 2.2 (Chen, et al. 5 ) 

Let {A*} denote the eigenvalues, and assume that 

|Afc| — oo. 


Then 


— c 

lim Re A* = - — 
k La 


( 8 ) 


Proof 

See Chen et a/. 5 for a proof. 

The authors of Chen et al. however do not appear to offer a proof of the fact that the eigenvalues {A*} are 
nonfinite in number. The fact that the resolvent is compact is not adequate to establish this; the compactness only 
assures that if nonfinite in number then {A*} can be arranged so that 

|Afc + i| > |A fc | 


and 


\X k \ — *> oo as k — ♦ oo. 


For proving the fact that eigenvalues are denumerably infinite we can indicate a general technique. 


Theorem 2.3 

The eigenvalues {A*} are denumerably infinite and such that 



< oo. 


( 9 ) 


Proof 

From (6) we see that for each a, the eigenvalues are the zeros of the function 

d( A) = (1 + Cosh 7 L Cos 7 L) - i 7 (Sinh 7 L Cos 7 L + Cosh 7 L Sin 7 I). (10) 

As a power series expansion will show, this is an entire function of the complex variable A. Moreover it is of 
exponential type, of order ±, and of completely regular growth. Further we can calculate that 

/i(0) = lim log |d(re' e )| = max(|Sinf|, |Cosf|). 

r— »oo 


5 



Let n(r) denote the number of zeros 
R.P. Boas (see Levin 7 ) we have: 


of <f(-) in the circle of radius r centered at zero. 

,im nTi = r m de > 0 

r — *oo r 1 '* 4 7T Jo 


Then by the theorem of 


Hence 

lim n(r) = oo, 

r— * oo 

or, the number of zeros is not finite. Moreover the function is of class A (see Levin 7 for the definition) since 

rR log|<*(s)d(-s)| 


sup 
R> 0 


J 


1 + s ' 2 


ds < Md < oo. 


The result (9) is a consequence. Q.E.D. 


Remark 

Applying Jensen’s Theorem we have 


J_ f 2 * d'{re' d ) 

2tt J 0 d(re i6 ) 


re ' 6 d 6 


n(r). 


We can acutally compute this as a quick means of locating eigenvalues. There is a jump of 2 corresponding to each 
eigenvalue and its conjugate. This is shown in Figure 1 for 



Leu 


Eigenfunctions 

The eigenfunction corresponding to the eigenvalue A* is given by 


= Afc 


<Pk 

^k<t>k 


<fi k (s) — c*:(Cosh 7 fc s - Cos 7 fc s) + dfc(Sinh7 fc s - Sin 7 k s) 


H{ A*) 


Cfc 

dfc 


= 0 


c fc = (Cosh7 fc L - Cos 7 fc L); d k = -(Sinh7 fc L - Sin7 fc L) 


where 
where 

or, we may take 
or, 

<t>k(s) = j4fc[(Cosh 7 fc L - Cos 7 fc L) (Cosh 7 fc s - Cos7 fc s) - (Sinh7 fc L - Sin7 fc i)(Sinh7 fc s 
Correspondingly: 

<t>' k (L) = 2A fc 7 fc (Cosh7 fc L - Cos7 fc L)Sin7 k L. 

The coefficient A k may be chosen for appropriate normalization. For example we may make 

ll*fcll= 1- 

Note that J k is an eigenvalue of A' and the corresponding eigenvector is: 

**(•) 

Vic = B k I 7 k {L) 

~^k<Pk(') 


- Sin7 fc s)]. (11) 
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where B k is again a “normalization” scalar. Note that 


[**, *fc]E 



= 4cckdkA k B k lk 

¥ o. 



7n & CoS 7n ^ 


Cosh7 fc L Sin7 fc L 


■f 


Ak&k 

Sinh7 n 5 Sin7 n 5) ds 


In particular we may choose A k> B k so that 

[*k,*k]E= I- ( 12 ) 

Further using a result of Gohberg and Krein 10 (we omit the details) we can establish that with the 

normalization (12) actually yield a R.iesz basis for £ . In terms of this basis we have the ( modal ) expansion for 
the solution of ( 2 ) 

oo 

Y(t) = ^a fc e A ‘ £ $ fc (13) 

1 

where 

at = (V'(O), 'I'*]*; 

and as an easy byproduct, using ( 8 ), we see that the semigroup generated by .4 is exponentially stable (established 
in Chen et ai by different arguments). 


Root Locus 

Let us consider how the eigenvalues behave as the gain a is varied. For this purpose it is convenient to define 

d( A;q) = M( A) -I- — AT (A) 


where 

M (A) = 1 + Cosh 7 L C 0 S 7 L 

N(\) = — — (Sinh 7 L Cos 7 L 4- Cosh 7 L Sin 7 L). 
cv 

Because of the analytic dependence of <f(A; a) on a, we can invoke the theory of algebraic or algebroidal functions 8,9 
and note that 

d(A(a);a) = 0 

will define A (a) as a multivalued analytic function of a with isolated singularities, if any. In particular this allows 
us to define the sequence {Ajt(a)}, k = 1,2, .. . such that 

A fc (0) = , Hk = (2k- 1)$ + e k 

L*V 


(the “clamped-free” beam modes) and 


lim Afc(a) 


i(kn — e'k) 


2 


l 2 v 


(“clamped-rolling” modes) and the real root 
is such that 


Ao(a) 


lim A 0 (q) = 0, lim A 0 (a) = - 00 . 

a— *00 a — >0 


A plot of the locus of the real root is shown in Figure 2. Moreover 


Afc(u) 


-1 


N(X) 


o' M'(A) + £JV'(A) 


A=A k (a) 
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In particular 

A' fc (0) = 

-1 N(A) I 
cv M'(A)| a=Aji(0) 


We can show that 

±(Re\ k (a)) 

da 

= $(£)■ “=° 

( 14 ) 



c 

= — — , a = +00 




2 a 2 L 



■j- (Im Afc(a)) > 0 . 
da 

A root locus of the first mode is shown in Figure 3 . The damping (= |ReA fc |) increases with the gain until a 
critical value of the gain is reached and thereafter decreases to zero. Note that by virtue of ( 14 ) we have actually 
“proportional damping” for small gain. A plot of the critical value of the gain versus the mode number is given in 
Figure 4 . 


Resolvent 

Let us now return to the resolvent — or solving ( 3 ). We note that 

0(A,s) = — — r f (Sinh7(s-<7) - Sin7(s - a)) v 1 (/i2(^) + ^i( a )) <*<7 

2 r Jo 

is a “particular” solution of 

a 2 u 2 f\ + /r = v\h 2 +\h x ) 

such that 

/i(0) = /i(0) = 0. 

Hence we can express the solution /i(A t s), where we have included A to indicate the dependence on A, as: 
= g(X, 5 ) + a(A)(Cosh 75 - Cos 7s) 4- 6(A)(Sinh7s - Sin7s), 0 < s < L 


where the coefficients a (A), 6(A) are determined from 


°(A) 

6(A) 


H( A)" 1 


q 6 - aXg'(X, L ) — cg"(X } L) 


where the primes again denote derivatives with respect to the variable s . Hence letting 


/i l i(A) 6. 12(A) 

H( A) = 

6^21 (A) 6 , 22 ( A ) 

and defining 

/122(A) — ^ 12 (A ) 

H(X) = 

-/121(A) * 11(A) 


so that 

H(X)H(X) = D{X)I = H{X)H{X) } 


a(A) — 1 , [^22(A)(q6 — aXg(X, L) - cg"{ X, L)) + * i2(A)p'" (A, L)\ 

D(\) 

5(A) = — l-/i2i(A)(a6 - aXg(X, L) - cg"(X,L)) - hn(X)g n, (X, L)] . 

D\ A) 
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We can cast the Green's function in the form: 


/ i(A.s) 


-i: 


W - ♦ Lj ^r 1 *“> * 


ab 

cW 


[/i 22 (A)(Cosh 7 « - C0S73) - /» 2 i(^)(Sinh 7 s - Sin7s)l 


K(X\( 7 ,s) = (Cosh 7 s - Cos 7 s) 


^2M _ ~~T h 22^j Cosh 7 (L - <r) + ^2 /i 12 + -^ 5 - * 22 ^ Cos -)(L - cr) 


2c/i22 


(Sinh 7 (L — <r) + Sin 7 (I -a)) 


4- (Sinh 75 — Sin 7 s) 


-^2a_U2i _ 2ftu ^ Cosh 7 (L- c) + ( ^f 21 - ~ 2 h„^Cos 7 (L-a) 


+ (Sinh 7 (L — a) 4 Sin 7 (L — a)) 

7 


3 < a 


(15) 


h = m(/i2 4- A/ii) 

£>(A) = — 2 A 2 i / 2 [ 07(1 4 Cosh 7 L C 0 S 7 L) 4- Aq (Sinh 7 L C 0 S 7 L 4 Cosh 7 L Sin 7 L)] 


Finally 


Note that setting a = 0 in (15) we get the Green’s function for the clamped/ free- free beam. In particular 



hi 


/ i(A,*) 

K(\,A) 

b 

= 

/{(A. 0) 


/l2 




n{o , .4) 

hi 

MO) 



Kh 2 

(Kh 2 )( 0) 

4 a/u(0) 

L=x 

c 

L 

c 


h .2 


-Ai 


0 


where Kh^ is the function given by 


— f (L — cr)/i2 (a) da 4- — f {L — a)/i2(cr) da, 0 < 5 < L. 
C Jo c Jl 
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Figure 1: n(|A|) Zeros. 
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Figure 2: Deadbeat Mode (Real Eigenvalue). 



Figure 3: Root Locus: First Mode A = a + iu>. 
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Figure 4: Critical Gain vs. Mode. 
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